Neutrinoless double /3-decay and neutrino mass hierarchies 
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> ■ I. INTRODUCTION 

o ■ 

, The discovery of neutrino oscillations in the atmos^eric Super-Kamiokande experiment ij , in the solar SNO 

' experiment yL in the reactor KamLAND experiment Q in the accelerator K2K experiment jj| and other neutrino 

] experiments la, 0, Q, '^'^s of the most compelling evidence in favor of new physics beyond the Standard Model. 

10 ' All existing neutrino oscillation data, with the exception of the data of the short baseline accelerator experiment 
LSND Q (the LSND result will be checked by the running MiniBooNE experiment jlOj soon) are described by the 

. three-neutrino mixing scheme 



Ph. I'lhix) ^'^Uii ViL{x); l = e,iJ,,T. (1) 

■ ' i—l 



Here, Viix) is the field of the neutrino with mass rrii {i — 1, 2, 3) and vil{x) is a flavor neutrino field which enters into 
the standard charged and neutral currents 



I I 
U is the unitary Pontecorvo-Maki-Nakagawa- Sakata(PMNS) [Hill 

mixing matrix. For massive Dirac neutrinos the 

PMNS matrix in the standard parameterization has the form 

C12C13 S12C13 5136"*'^ 

= I -S12C23 - ci2S23Si3e*'' C12C23 - si2S23Si3e''' S23C13 | • (3) 

S12S23 - Cl2C23Sl3e** -C12S23 - 5120235136"^ C23C13 

Here Sy — sin&y , Cy — cos^y , 0ij (i < j) is the neutrino mixing angle and 5 is the CP violating phase. 
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From the analysis of the Super-Kamiokande atmospheric neutrino data for the neutrino mass squared difference 
Am^g and the parameter sin^ 2^23 it was obtained (l| : 

best fit : Am^3 = 2.1 10"^ eV^ sin^ 26I23 = 1.00, 



''23 

90% C.L. : 1.5 10-3 eV^ < Ani^g < 3.4 10"^ eV^ sin^ 26*23 > 0.92. (4) 



The global analysis of the data of the solar neutrino experiments and KamLAND experiment yields the following best 
fit values and 90% C.L. ranges of the relevant neutrino oscillation parameters 0: 

best fit : Am^2 = 7.9 10"^ eV^ tan26li2 = 0.40, 

90% C.L. : 7.4 10"^ eV^ < Am^2 < 8.5 10"^ eV^ 0.33 < tan^ 6(12 < 0.50. (5) 

Notice that neutrino mass-squared difference is determined as Am^^ = to| — mf . For the angle ^13 only upper bound 
is known. From the exclusion plot obtained from the data of the reactor experiment CHOOZ 0,^3 '^^ have 

sin^ 6*13 < 5 10^2 (90% C.L.). (6) 

The CP-violating phase S remains undetermined. A recent global analysis of the oscillation data lead to the following 
bound: sin^ 6I13 < 0.9lo.glO-2 (95% C.L.) 

At present the structure of the neutrino mass spectrum is not known as well. Two types of spectra are possible: 

1. Normal spectrum: 

rrii < 7712 < TO3; Am^2 ^ Am|3. (7) 

2. Inverted spectrum: 

m3 < mi < TO2; Ami2 ^ I^™i3l (8) 

We note that it is common to label neutrino masses differently in the case of the normal and the inverted spectra. 
For both spectra we have m2 > mi. But in the case of the normal spectrum m^ is the mass of the heaviest neutrino 
and in the case of the inverted hierarchy 777,3 is the mass of the lightest neutrino. This convention allows to keep the 
same notation of the mixing angles for both spectra. Existing oscillation data are compatible both with normal and 
the inverted spectra. 

The lightest neutrino mass mg = 7711(7713), which determines the absolute values of neutrino masses, is currently 
also unknown. From an analysis of the data of the Mainz 16| and Troitsk |17| tritium experiments it was found 

Too < 2.3 eV. (9) 

A more stringent bound on the sum of neutrino masses can be found from the measurement of the matter power 
spectrum P{k). Depending on the data which were taken into account, the cosmological upper bound on the sum of 
neutrino masses was obtained as (see 0, 0| and references therein) 

^TO, < (0.5 - 1.7) eV. (10) 

i 

An important evidence that masses and mixing of neutrinos are of a nature beyond the Standard model (SM) would 
be that massive neutrinos are Majorana particles. If i^i are Majorana particles 

1. Neutrino fields i'i{x) satisfy the Majorana conditions 

,y^{x)=,y,{x), (11) 
where i'f(x) = Cvf{x) is the conjugated field (C is the charge conjugation matrix). 

2. The neutrino mixing matrix has the form |20| 

U = U'^S{a) (12) 

where S{a) is a diagonal phase matrix. In the case of three neutrino mixing the matrix S{a) is characterized 
by two Majorana CP-violating phases. The matrix S{a) can be presented in the form 

S^k ^ e"'^ 6,k; 03 = 0. (13) 

The unitary matrix C/^, which is characterized by the three mixing angles 6'i2, ^23, ^13 and one phase 5, was 
already introduced in Eq. 
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If in the lepton sector CP invariance holds, for the Majorana mixing matrix we have |21| 

Uh^U;,!!,, (14) 
where iji — ±i is the CP parity of the Majorana neutrino Vi. The condition (|14|l can be presented in the form 

Ul,^\Uu\^e''-^/^^P% (15) 

where pi = ±1. 

Investigations of neutrino oscillations in vacuum and in matter do not allow to distinguish massive Dirac from 
massive Majorana neutrinos [20l 123 . |2^ . In order to reveal the Majorana nature of Vi it is necessary to study 
processes in which the total lepton number is violated. Because the standard electroweak interaction conserves helicity 
the probabilities of such processes are proportional to the squares of the neutrino masses, and, consequently, they are 
strongly suppressed. The best sensitivity on small Majorana neutrino masses can be reached in the investigation of 
neutrinoless double /3-decay (0z//3/3) of some even-even nuclei. 



II. NEUTRINOLESS DOUBLE /3-DECAY 



In the case of Majorana neutrino mixing and the standard electroweak CC interaction the Oi'/J/J-decay, 

{A, Z) ^ {A,Z + 2) + e- +e-, (16) 

is a second order process in the Fermi constant Gp with virtual neutrinos. The half- life of the process is given by the 
following general expression .24ij 

|m^^P|MO'^(A,Z)pG°^(£;o,Z). (17) 



-1/2 

Here 



'^Pfi = X! ^ez mi (18) 



is the effective Majorana mass, M^'^{A,Z) is the nuclear matrix element (NME) and G^^ {Eq, Z) is a known phase- 
space factor [Eq is the energy release). Let us stress that the NME is determined only by nuclear properties (its 
dependence on the small neutrino masses can be safely neglected). 

After the discovery of neutrino oscillations, the search for neutrinoless double beta decay (Oi^/3/3-decay) became one 
of the most fundamental problems of neutrino physics. Observation of this process would be the proof that massive 
neutrinos Vi are Majorana particles l^^. Furthermore the observation of Qvj3(3-Aec&y will allow to reveal the type of 
the neutrino mass spectrum, to determine the mass of the lightest neutrino and, possibly, Majorana CP phases. 

The most stringent lower bounds on the half- life of Oi//3/3-decay were obtained in the Heidelberg- Moscow ^^Ge 
and CUORICINO "Oxe |23 experiments: 

T{'/2(™Ge) > 1.9 • 10^5 years, T'^j2(^^"Te) > 1.8 ■ lO^^ years. (19) 



Using recently calculated nuclear matrix elements with significantly reduced theoretical uncertainties [28|, |29j from 
these data the following upper bounds for the effective Majorana mass can be inferred 

|"T-/3/3| < 0.55 eV (Heidelberg — Moscow) 

\mf3p\ < 1.1 eV (CUORICINO). (20) 

The Heidelberg group, which includes a few authors of the Heidelberg-Moscow collaboration, recently claimed |30| 
evidence for the Oi^/3/3-decay of ^^Ge with T^^^ = (0.69 - 4.18) • 10^^ years at the 4.2cr confidence level. Using the NME 
obtained in Ref. [28, 29], from this data one finds for the effective Majorana mass the range 0.37 eV < Im^p] < 0.91 eV. 
The claim made in 30J was re-examined and critized by different authors and in particular by the Moscow 
participants of the Heidelberg-Moscow collaboration [s^l ■ The GERDA I experiment , now at preparation at Gran 
Sasso, will be able to check relatively soon the claim made in |30| . 

The effective Majorana mass given by Ea. H18|l is determined by the values of the neutrino masses rrii, which for the 
case of the normal neutrino mass spectrum are given by 



m2 = y ml + Aml2, ma = ymf -I- Amfj + Am^g, (21) 
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and by the matrix elements U^j^ , which in the standard parameterization take the form 

C/^i = cos^ cos^ 012 e^*"^ , [7^2 = cos^ sin^ 9^2 e''"'^ , C/'a = sin^ ^13 e^'"^ , (22) 

where are Majorana phases. 

The values of the neutrino masses depend on the lightest neutrino mass mo = 7711(7713), on the neutrino mass 
spectrum and the neutrino mass squared differences Ami2 and Am^g ( jAmigl), which are known from neutrino 
oscillation data (see and |(SJ) ) . The value of the parameter sin^ 0i2 is deduced from the analysis of the solar and 
KamLAND data (see (O). The lightest neutrino mass toq and the CP Majorana phases are unknown and will be 
considered as free parameters. 

In the case of the normal mass hierarchy, 

TOl ^ 7712 <C "73, (23) 

the lowest two neutrino masses and the effective Majorana mass |777^^| have the minimal values. If we neglect the 
contribution of 7711 to H18|) . for the upper bound of |777^^| we get 



|to/3/3| < (sin^ 6*12 J A777^2 + sin^ ^i3 V A777^3). (24) 



The contribution of the first term to \mi3/3\ is small because of the smallness of -y/ Aml2- The contribution proportional 
to the "large" ^ ATn^g is suppressed by the smallness of the parameter sin^ ^13. With (0), jnj and (|24(l for the upper 
bound of the effective Majorana mass we find 

\m(3(3\<6A- 10-^. (25) 

In the case of the inverted hierarchy 

7713 < mi < (26) 

we can safely neglect the contribution of the lightest mass m3 to the effective Majorana mass. For |m^/3| we have the 
following expression 



|m;3;3| ~ ^lAm^gl (l-sin^2 0i2 sin^a2i)'/', (27) 
where a2i = a2 — ai. The only unknown parameter in (|27() is sin^ 0:12. From (|27|l we have 



Amfgl cos26'i2 < |m/30| < ^lAmfgl. (28) 

The bounds in 1)281) correspond to the case of the CP invariance in the lepton sector (the upper bound corresponds 
to pi = p2 and the lower bound corresponds to pi = —p2-) From and (|28|l we find that in the case of the 
inverted hierarchy the value of Im^j^j must lie in the range 

I.O-IQ-^ <\mf30\< 5.5-10-^ eY, (29) 

In the case of the quasi-degenerate spectrum of neutrino masses, 



mi < m2 < mg; m^ ~ mo 3> \/ Am23, (30) 



the effective Majorana mass is given by Eq. (|27|l in which the replacement -^/jAmigl — > 7710 must be performed. 

Thus, in the case of the quasi-degenerate spectrum the effective Majorana mass depends on two parameters: nio 
and sin^ Q!2i. For the common neutrino mass we have 

cos Z (712 

The current bound on mo, obtained from the measurement of the high energy part of the /3-spectrum of tritium, is 
given in 0. The future KATRIN tritium experiment will be sensitive to a value of mo ~ 2 • 10~^ eV. 

There are many models of neutrino masses (see e.g. Refs. |35l l36l ISTj ). The see-saw mechanism of neutrino mass 
generation j38l | is considered to be the most plausible one. This mechanism is based on the assumption that violation 
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TABLE L Sensitivities of future Oi/13/3-deca.y experiments to the effective Majorana neutrino mass |m;3^| calculated with the 
RQRPA nuclear matrix elements M^'''^{A, Z) of Ref. |2fl|]. For the axial coupling constant gA the value gA = 1-25 was assumed. 
T^J^"^^^ is the maximal half-life, which can be reached in the experiment and |m^/3| is the corresponding upper limit of the 
effective Majorana neutrino mass. 



Nucleus 


Experiment 


Source 


-'1/2 


-exp r 1 

[yr] 


Ref. 


M°''{A,Z) 


Imp/jl [eV] 


^«Ge 


GERDA(I) 


15 kg of '^"'■Ge 


3 


^q25 




2.40 


0.44 




GERDA(II) 


100 kg of ""■Ge 


2 


^(,26 


133] 


2.40 


0.17 




Majorana 


0.5 t of ™ Ge 


4 


10^^ 


M 


2.40 


0.038 




SuperNEMO 


100 kg of ""''Se 


2 




m 


2.12 


0.091 




MOON 


3.4 t of ""*A/o 


1 




170] 


1.16 


0.058 




CAMEO 


1 t of CdW04 crystals 




^(,26 


no] 


1.43 


0.14 


"Ore 


CUORE 


750 kg of TeOa 




10^^ 


[73] 


1.47 


0.047 




XMASS 


10 t of liq. Xe 


3 


10^6 


f70] 


0.98 


0.12 














0.73 


0.17 




EXO 


1 t '="''Xe 


2 


10^7 


[75] 


0.98 


0.048 














0.73 


0.064 



of the total lepton number is at a large scale and connects the smallness of the Majorana neutrino masses with heavy 
right-handed Majorana particles. The existence of such particles provide a natural framework for the explanation of 
the baryon asymmetry of the Universe (see Ref. H^)- 

In the framework of the see-saw mechanism the degenerate neutrino mass spectrum requires a fine-tuning which 
includes the Dirac mass matrix and the right-handed Majorana mass matrix. 

The inverted hierarchy of neutrino masses requires a specific lepton symmetry like a global gauge symmetry which 
provide conservation of — L^^ — Lr- In order to explain existing data this symmetry must be broken. In such a 
framework it is difhcult to reconcile the large 623 with the not so maximal mixing angle 0i2 and the small mixing angle 
^13 (35I [s^ . Neutrino mass hierarchy is a natural spectrum in the case of the see-saw mechanism. Such a spectrum is 
realized in the case of SO(IO) ^3 another GUT models which connect quark and lepton sectors. The problems 
of the CP phases and the renormalization group effects in the SO(IO) GUT were discussed in 0. 

In this paper we calculate the half-lives of the Oi^/3/3-decay in the case of the normal and inverted hierarchies of 
neutrino masses. Neutrino masses for such spectra are very small. For the case of the normal hierarchy we have 

m2 ~ 9 • 10"^ eV; mg ~ 5 • lO'^ eV; mi < 1R2 (32) 

For the case of the inverted hierarchy we obtain 

mi ~ 1712 — 5 • 10~^ eV; ma <^ nii. (33) 

Effects of such small neutrino masses cannot be observed in tritium and other /?-decay experiments in a foreseeable 
future. Future data on the distribution of clusters of gala xies and gravitational lensing data, however, will be sensitive 
to the following value of the sum of neutrino masses IT9I li^ 

^m, ~ 3-10-2eV. (34) 

i 

Thus, apparently, future cosmological measurements can probe neutrino mass hierarchies. 

Several future experiments on the search for Oj//3/?-decay will be sensitive to values of the effective Majorana mass 
in the range (|29|l .which corresponds to the inverted mass hierarchy (see Table IJ). 

Taking into account the theoretical plausibility of the normal hierarchy of neutrino masses we believe that it 
is worthwhile to consider the 0uPf3-decay in detail for this type of neutrino mass spectrum. Recently, important 
progress in the evaluation of the Oi^/J/J-decay NME's for the ground state transitions of nuclei of experimental interest 
was achieved 28, 29]. We shall use these new results to calculate expected half-lives of Oi//3/3-decay for both neutrino 
mass hierarchies. 
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III. NUCLEAR MATRIX ELEMENTS 



From the measurement of the half-Hfe of the Oi^/3/3-decay only the product |m^^| |Af°'^(A, Z)\ can be determined. 
Thus, without accurate calculation of nuclear matrix elements, it is n ot p ossible to reach qualitative conclusions about 
neutrino masses and the type of neutrino mass spectrum |4a. lil. 1451 149. 1471 148| . 

The calculation of the Oi^/3/9-decay matrix elements is a difficult problem because ground and many excited states 
of open-shell nuclei with complicated nuclear structure have to be considered. In the calculation of the Oi^/3/3-decay 
NME's the nuclear shell model (NSM) and the proton-neutron quasiparticle random phase approximation (pn-QRPA) 
or extensions to it 49] are used. 

These two approaches are significantly different. The NSM is limited to a set of single-particle states in the vicinity 
of the Fermi level. Thus, configurations with only small excitations are considered. These excitations are correlated 
in all possible ways. The open problem is the effect of single-particle states further away from the Fermi level, which 
is neglected. The large NSM spaces face the problem of diagonalization of large matrices and the construction of good 
effective interaction. 

From the available shell-model calculations of the NME of the Oi//3/3-decay the most advanced are those of the 
Strasbourg group [s^l, which appeared about ten years ago. From that time, in spite of significant progress in 
computer speed and memory no large-scale NSM calculations of the Oi//3/3-decay NME have been published. 

The pn-QRPA method takes into account many single-particle states including states relatively far from the Fermi 
surface but with correlations which are of the specific simple type. Many extensions of the pn-QRPA have been 
proposed with the aim to improve the many-body approximation scheme. The so-called renormalized pn-QRPA 
(pn-RQRPA) scheme takes into account the Pauli exclusion principle by improving the quasi-boson approximation 
|5l|. In this method for the evaluation of the commutators of bifermion operators the QRPA ground state is used 
(instead of the uncorrelated BCS ground state). This refinement of QRPA approach allows to reduce significantly 
the sensitivity of the results to the renormalization of the particle-particle interaction of the nuclear Hamiltonian. 
The fact that pn-RQRPA is an improvement of the pn-QRPA approach has been confirmed by the studies performed 
within a schematic model js^ . 

The pn-QRPA and its extensions remain a popular technique of the calculation of Oz^/3/?-decay NME's. However, 
various implementations of the QRPA introduced by different authors have produced a spread of results with a 
factor of 2-3 and even more difference in the NME '53\. Some authors simplified this problem by assuming that the 
published range of calculated NME's defines a plausible approximation to the uncertainty in our knowledge of the 
matrix elements 113, 01 • We do not share this position. We believe that it is not appropriate to consider all calculated 
0://3/?-decay NME's at the same level. The correct procedure is to understand the difference among various QRPA-like 
calculations and the origin of contradictions between different results. In Ref. p9l | a list of main reasons leading to a 
spread of the pn-QRPA and the pn-RQRPA nuclear matrix elements was presented. 

One of the most important factors of the QRPA calculation of the NME's is the way how the particle-particle 
strength of the nuclear Hamiltonian gpp is fixed. When the early QRPA calculations were performed only a limited 
information on half-lives of the 2i^/3/3-decay was available. At that time gpp was fitted to existing data on single P~- 
transitions or alternatively gpp was chosen to be equal to unity. Nowadays, the half-lives of 2i'/3/3-decay of ten nuclei 
has been measured. Recently, it has been shown that by adjusting gpp to the 2^/3/3-deca y ra tes we can significantly 
eliminate uncertainties associated with variations in QRPA calculations of decay rates |28l l29j| . In particular, the 
results obtained in this way are essentially independent of the size of the basis, the form of different realistic nucleon- 
nucleon potentials, or on whether QRPA or RQRPA is used. Furthermore, the matrix elements are shown to be also 
rather stable with respect to the possible quenching of the axial vector constant gA- 

The procedure proposed in ,2Si, critically analyzed in [s^. The author's conclusion was that fitting of 

gpp to (3^ (or electron capture (EC)) and single /?~-decay of the ground state of the intermediate nucleus is a more 
meaningful procedure. This criticism has been refuted in Ref. ^9|. In particular, it was shown that there is no reason 
to give preference to the lowest state of the intermediate nucleus. The (3~ and / EC matrix elements move with 
gpp in opposite directions, what makes it difficult to adjust gpp by choosing one of them. It is preferable to use sum 
of the products of amplitudes, i.e., the 2j//3/3-decay half-life. It was also noticed that practically for all multipolarities 
significant amount of strength is concentrated up to 10-15 MeV and that the contributions of the 1+ multipole to the 
2v(3(3- and Oz^/?/3-decay matrix elements are correlated. Thus, there is no reason to choose any one particular state or 
transition for adjustment. It is worth also mentioning that only for three double /?-decay nuclear systems of interest 
(A — 100, 116, 128) the ground state of the intermediate nucleus is just the 1+ state. Thus, the procedure of how to 
fix gpp proposed in Ref. ^5| is not only disfavored but also strongly limited. These and other arguments of Ref. 
clearly favor the procedure of the 2v (3 j3-decscy fitting rather than the procedure of the beta-decay fitting. 

A discussion concerning the previous QRPA and RQRPA calculations of the 0^/9/3-decay NME's, in which different 
assumptions and approximations were made, can be found in p9| . It was shown that in most, albeit not all, cases 
the differences among them can be understood. Attention was also paid to the fact that the results of (5^ differ 
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significantly from those of |23,l29!| in spite of the fact that the same procedure of adjustment of gpp was used (for the 
most important 1+ channel). Contrary to the NME of [2^ the NME calculated in [SQl] strongly depend on the 
size of the model space. In particular the levels lying far from the Fermi surface severely influence the decay rate. 
Actually, there is no explanation of this fact and nobody else reported such a strong effect. 

The calculations performed in Ref. can be considered as the most reliable QRPA and RQRPA calculations 
of the NME of the Oz//3/3-decay. We stress that the NME calculations in which the dependence of the results on the 
particular choice made is not discussed could not be considered on the same footing as those where these points are 
carefully explained (see e.g. Ref. [29|'). 

There are no doubts that further progress in the calculation of the Oi//3/3-decay NME's is needed. The nuclei 
of experimental interest for the investigation of Oi'/J/^-decay have been extensively studied by assuming spherical 
symmetry. However, many of the nuclei undergoing double beta decay are deformed {57| and it is important to study 
the effect of deformation on the Oz^/3/3-decay. Recently a new suppression mechanism of the 2i'/3/?-decay matrix elements 
based on the relative deformation of the initial and final nuclei has been found ||53]. The effect of deformation is large, 
if there is a significant difference in the deformations of the parent and daughter nuclei. The effect of deformations 
could also be large in the case of the Oi^/3/3-decay matrix elements. 

For further progress in the field it is also important to apply the methods of the calculation of the Oi^/3/3-decay 
matrix elements to the calculation of the matrix elements of related processes like charge-exchange reactions 59j, 
muon capture 60] and charged current (anti)ncutrino- nucleus reactions '6l'|. The observation of these processes 
might probe forbidden transitions which contribute to the Oi^/3/3-decay half-lives. Thus, it is essential to understand 
how the methods applied for the calculation of the NME of the Oz//3/3-decay can reproduce the observables related to 
these processes. For the improvement of many-body approaches it is also important to test them in the cases of exact 
solutions of solvable models which are as realistic as possible. 

The improvement of the calculation of the nuclear matrix elements is a very important and challenging problem. 
The nuclear matrix elements of the Oz^/3/3-decay cannot be related exactly to other observables. A complementary 
experimental information from related processes is highly required, but it cannot fully solve the problem of the 
uncertainties in the 0^'/3/3-decay matrix elements. This problem might be solved only by observation of the 0i>(3P- 
decay of at least three different nuclei ^ . This would be a model independent test of the theoretically calculated 
NME's. 



IV. THE EXPECTED Oi^/3/3-DECAY HALF-LIVES 

Here, we present the results of the calculation of the half-lives of the Oi//3/3-decay of following nuclei ^^Ge, ^^Se, 
iooMo, ^^®Cd, ^•^"Te and ^^^Xe many of which presumably will be investigated in future experiments. We will use 
nuclear matrix elements M^" {A, Z) obtained in Ref. [2^ . 

First discuss the case of the normal hierarchy of neutrino masses. Taking into account all terms for {mppl we have 
the following expression 

\mp(3\ = I cos^ 6*12 mi + sin^ 6*12 sj Amj^ e^'°"^ + sin^ 6*13 ^Amig e'^'"''\. (35) 

There are four unknown parameters in (|35|) . namely ™i,sin2 and two CP phase differences ai2 and ai^. First we 
shall assume that the lightest mass mi is very small and we can neglect its contribution to |to/3/3|. In this case the 
effective Majorana mass m^^ depends on two free parameters: sin^ ^13 and 023. 

The current 90 % CL CHOOZ bound on the parameter sin^ ^13 is given in ©. A significant improvement on the 
value of this parameter is planned to be achieved in future experiments. The expected sensitivity of the Double- 
CHOOZ reactor experiment "bS] is sin^ ^13 ~ 7.5 ■ 10~^. In the accelerator T2K z/^ experiment [6^ a sensitivity 

of sin^ 6'i3 ~ 1,5 • 10""^ is planned to be achieved. Thus, we will consider for the parameter sin^ ^13 the range 
(10-3 - 5 • 10-2). 

In Fig. [T]we present the half-life Tl'^^^A, Z) as a function of sin^ 6I13 for the nuclei ^^Ge, ^^Se, ^°°Mo, "^Cd, ^^°Te 

and -'^^^Xe. The regions with solid line boundaries (dashed line boundaries) were calculated with the best fit values 
( 90% C.L. values) of the neutrino oscillation parameters Am23, Ami2 and sin^ 6*12 (see Eqs. Q and (O). The 
boundaries of these regions correspond to the case of CP-conservation. The lower (upper) boundaries correspond to 
cos2q;23 = 1 (cos2a23 = —1), i-e., to the case of same (opposite) CP parities of the second and third neutrino. For 
the case of the Oz^/3/3-decay of ^^^Xe the allowed regions in Fig^are larger than for other nuclei. This is connected to 
the fact that the corresponding NME is constrained to a range of values since the 2i//3/3-deca.y of this isotope has not 
been measured yet |^] . 

From Fig. ^it follows that for sin^ 6*13 < 10"^ the calculated half-lives practically do not depend on sin^ 6*13. For 
such small values of sin^ ^13 the dominant contribution to the effective Majorana mass m^p comes from the solar term 
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FIG. 1: The neutrinoless double beta decay half-life T^'j,^ for nuclei of experimental interest as function of the parameter 
sin^ ^13 . The case of the normal hierarchy of neutrino masses is considered and the lightest neutrino mass is assumed to be 
negligibly small. The region with solid line (dashed line) boundaries corresponds to the best fit (90% C.L.)of the neutrino 
oscillation parameters U Q. The calculations were performed by using recently evaluated nuclear matrix elements with 
significantly reduced theoretical uncertainty |29||. The vertical line indicates the current upper limit on sin^ &13 set by the 
CHOOZ experiment. 



(the second term on the right hand side of Eq.(|23J)). From Fig. Qwe also deduce that if the limit on sin^ ^13 will be 
pushed further down by the Double-CHOOZ and T2K experiments the expected half-lives of the Oi^/?/3-decay of ®^Se, 
iooMo, "^cd, i30Te and ^^Sxe will be about lO^^ - 10^° years. 

In the Table Hll we present the minimal values of T^^^i^A^Z) in the case of the neutrino mass hierarchy with a 

negligibly small lightest mass. They are given for three possible values of the parameter sin^6'i3: 0.05, 0.01, 0.001. 
These values are expected to be slightly increased if the limit on the parameter sin^ 6113 will further decrease (see 
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TABLE II: Normal hierarchy of neutrino masses: The neutrinoless double beta decay half-life T°J^{A,Z) of '^''Ge, *^Se, ^^Zr, 
^°°Mo, ""^Cd, ^^*T^ ^^°Te ^^'^Xe and ^^"Nd The results are presented for three values of 613 from the allowed range sin^Ois < 
0.05 of CHOOZ |l4| . The best fit and 90% C.L. values of T^^2 were calculated by assuming equal values of Majorana CP- 
violating phases, i.e., the lowest allowed values for T^^jC^)-^) ^-rs given. The 0://3/3-decay nuclear matrix elements of Ref. |29| 
are used. 

Normal hierarchy: T°';'2{A,Z) [years] 
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Fig. ^ . These values can be confronted with the sensitivities on half-hves of future 0^/3/?-decay experiments given in 
Table HJ From the comparison of the calculated half-lives with the experimental sensitivities we conclude that further 
experimental efforts will be required to reach lower limits comparable to the predicted values in the case of the normal 
hierarchy of neutrino masses. 

Up to now we neglected the contribution of the lightest neutrino mass mp = toi to |to^/3|. From %VZ\ it follows that 
TO2 ~ 0.2 ma. If we (arbitrarily) assume that mi ~ 0.2 TO2 it follows that the modulus of the first term in (j35|l is about 
half of the modulus of the second solar term. Thus, the contribution of the mi-term to |m/3^| can be sizable. In Fig. 
|21we present three-dimensional plots for the half-life of the Oj^/3/3-decay of ^^Ge as a function of mi and sin^ ^13 under 
different assumptions for the relative CP Majorana phases. In the case of other isotopes similar results are expected. 
From Fig. [3 it follows that for a small values of sin^ ^13 < 0.01 the pronounced hill region with relatively large values 
of the half-lives is not more accessible. 

In the case of the inverted hierarchy of neutrino masses the upper and lower bounds on the effective Majorana mass 
\mpj3\ are well determined by the known oscillation parameters (see (I28II '). In Table HTll we give the corresponding 
ranges for the half-lives of relevant nuclei calculated with the nuclear matrix elements of "2^ by taking into account 
the best fit and 90% C.L. values of the parameters sin^ 612 and jAm^gj. In the same Table the predicted ranges of the 
half-lives are compared to the current experimental lower limits on the Oi^/3/3-decay half-life from experiments with 
saturated sensitivity and from two running experiments (NEMO 3 [6^ and Cuoricino |23|) together with half- lives 
sensitivities of proposed future experiments. We note that the designed future Oi^/3/3-decay experiments allow different 
strategies. Some of them plan to proceed with smaller steps forward and some of them prefer to make large steps. The 
time-scale for these experiments except the GERDA I |33| , which main task is to confirm or rule out the recent claim 
of evidence for the Oi^/3/3-decay of ''^Ge [s^l, is not determined yet. From Table UllI we conclude that the sensitivities of 
the next generation QvfSP-Aecay experiments will apparently allow to probe the inverted hierarchy of neutrino masses. 
Finally, in Fig|31we present the expected life-times of the Qv (3 (3-(iecay of ^^Ge, ^^°Te, ^^^Xe and other nuclei calculated 
for 90% CL allowed values of the parameters. The dependence of the 0:^/3/3-decay half-life on sin^ 612 is outlined. 
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FIG. 2: The neutrinoless double beta half-life T^^j^ of ™Ge as a function of the parameter sin^ Qn and the lightest neutrino 
mass mi. The normal hierarchy of neutrino masses is assumed and the best fit values of neutrino oscillation parameters are 
taken into account. The relative CP Majorana phases of neutrinos are denoted by 021 and aai. 



V. CONCLUSION 

The effective Majorana mass, which determines the half-life of neutrinoless double /3-decay, crucially depends on the 
character of the neutrino mass spectrum. All possible physical neutrino mass spectra (hierarchy, inverted hierarchy 
and quasi degenerate) are at present viable. However, in the framework of the see-saw mechanism, which is a plausible 
explanation of the smallness of neutrino masses, hierarchy is a favorable spectrum (see, |35l l36j|). Generically, the 
neutrino mass hierarchy naturally appears in GUT models like SO(IO) which unify quarks, leptons and neutrinos. 

Having in mind that the hierarchy of neutrino masses is a plausible neutrino mass spectrum, in this case we 
performed detailed calculations of half-lives of the Oi^/3/3-decay of several nuclei of experimental interest. We used 
nuclear matrix elements which were obtained in the recent most updated and most comprehensive RQRPA calculations 
[2^. We studied the dependence of the half-lives of 0i//3 /3-decay on sin^ ^13 and the lightest neutrino mass m\. The 
current ^3 hmit and future Ig^ B3| sensitivities to the value of the parameter sin^ 6'i3 were considered. The calculated 
lower limits on the half-lives of nuclei considered are listed in Table Hll As it is seen they are in the range of 10^* — 10^° 
years. The expected half-life sensitivities of the next generation of the Oj^/3/3-decay experiments are significantly lower. 

The future Oi^/3/3-decay experiments will probe the inverted hierarchy of neutrino masses which requires a symmetry 
in the neutrino mass matrix. Using the updated RQRPA calculations of NME [2^ we have calculated half-lives of 
Oi//3/?-decay of several nuclei in this case. 
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TABLE III: Inverted hierarchy of neutrino masses: The neutrinoless double beta decay half-life T^j2{^y^) calculated in 
RQRPA for "«Ge, ^^Se, ««Zr, i»«Mo, "«Cd, ^^^Te, "OTe ^^^Xe and ^^ONd 

C.L. ranges of the corresponding parameters. The Oi//3/3-decay nuclear matrix elements of Ref. 



Tip was calculated for best fit values and 90% 
1231 are used. The current 



experimental lower limits on the Ou/SfS-deca.y half-lives {T^J^ '^^^) are from experiments with saturated sensitivity except those 
denoted by the symbol *, which indicate experiments still taking data. H-M means Heidelberg-Moscow. 
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3: The neutrinoless double beta decay half-life 



for nuclei of experimental interest as function of the parameter 
9i2 ■ The case of the inverted hierarchy of neutrino masses is assumed. Conventions are the same as in Fig. QThe vertical 
dashed lines correspond to the boundaries set by the 90% C.L. values of 612 (3]. These results are not sensitive to ma for values 
below 10-2 el/ |4al. 
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